Ghost Imaging: What is quantum, what is not 
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We provide a unified treatment of ciassical and quantum Gaussian-state sources that unambigu- 
ously identifies which features of ghost imaging are strictly quantum mechanical. We show that 
ghost-image formation is fundamentally classical, with the image being expressible in terms of the 
phase-insensitive and phase-sensitive cross correlations between the detected fields. We then con- 
sider ghost-imaging scenarios with either phase-insensitive or phase-sensitive sources, where the 
former are always classical but the latter may be classical or quantum mechanical. We show that if 
their auto-correlations are identical, then a quantum source provides resolution improvement in its 
near-field and field-of-view improvement in its far field when compared to a classical source. 

PACS numbers: 42.30.Va, 42.50.Ar, 42.50.Dv 



Ghost imaging is the acquisition of an object's trans- 
mittance pattern by means of intensity correlation mea- 
surements. Its first demonstration utilized a biphoton 
source, thus the image was interpreted as a quantum phe- 
nomenon owing to the entanglement of the source pho- 
tons However, subsequent experimental @, H, 0] and 
theoretical [f], [f| considerations have demonstrated that 
ghost imaging can be performed with thermalized laser 
light. Whereas the biphoton requires a quantum descrip- 
tion for its photodetection statistics, thermal light does 
not, i.e., it can be regarded as a classical electromagnetic 
wave whose photodetection statistics can be treated via 
shot-noise theory. This disparity has sparked interest in 
determining whether the fundamental physics of ghost 
imaging require a non-classical source 

Bum nil- in 

this paper, we shall treat ghost imaging with classical 
and non-classical Gaussian-state sources to quantify the 
classical/quantum boundary. We begin by deriving the 
image — in terms of the detected fields' cross-correlation 
functions — for arbitrary two-mode, coherence separable 
Gaussian-state sources. We then compare image reso- 
lution and field-of-view for near-field and far-field ghost 
imaging using three types of sources: thermal light, clas- 
sical phase-sensitive light and quantum phase-sensitive 
light. We conclude with a discussion of the relevant 
physics in ghost imaging, emphasizing which aspects are 
classical and which are quantum. 

Consider the ghost imaging setup shown in Fig. [TJ 
Here, E s (p,t)e~ iuJ(,t and E R (p,t)e~ iuJot are a pair of 
scalar, positive frequency, z-propagating field operators 
normalized to photon-units, where ojq is their common 
center frequency, p is the transverse coordinate with re- 
spect to each one's optical axis, and 



[E m (Pi,ti), El(p 2 , t 2 )] = S(p 1 - p 2 )S(t 1 - t 2 



(1) 



for m £ {S, R}, are the non-zero commutators. Both 
beams undergo quasimonochromatic paraxial diffraction 
over an i-m-long free-space path, yielding measurement- 
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FIG. 1: (Color online) A simple ghost imaging setup. 



plane field operators 



E e (p,t) = J dp' E m (p',t-L/c 



k oe ik \p-p'\ 2 /2L 



i2nL 



(2) 



where (£, m) — {(1, S), (2, R)}, c is the speed of light 
and ko — u>o/c. The first field, Ex(p,t), illuminates a 
quantum-limited pinhole photodetector located at trans- 
verse coordinate p l . The second field, E 2 (p,t), illumi- 
nates an amplitude-transmission mask T(p), located im- 
mediately in front of a quantum-limited bucket photode- 
tector with sensitive region p £ A%. The product of the 
photocurrents from these sub-unity quantum efficiency, 
finite-bandwidth detectors is time averaged to estimate 
their ensemble-average cross correlation, C(p 1 ). This 
process is repeated, as p 1 is scanned over the plane, to 
obtain the ghost image of the object's intensity transmis- 
sion |T(p)| 2 . 

Without appreciable loss of generality, we shall assume 
that Eg and En are in a zero-mean, space-time coherence 
separable, jointly-Gaussian quantum state. Such a state 
is completely characterized by the field operators' phase- 
insensitive (normally-ordered) and phase-sensitive auto- 
and cross correlations. We allow the source fields to have 
both phase-insensitive and phase-sensitive cross correla- 
tions, but for simplicity we shall assume they carry no 
phase-sensitive auto-correlations. Because © is linear, 
Ei and E 2 are also in a zero-mean, coherence separable, 
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jointly-Gaussian state, and hence fully determined by 

(El( Pl ,h)E e (P2,t2)) =K^( Pl ,p 2 )RS e (t 2 -h) (3) 
{E 1 {p 1 ,t 1 )E 2 {p 2 ,t 2 )) = K { 1 p 2 \p 1 ,p 2 )R { $(t 2 -t 1 ), (4) 

for m,£ € {1,2}. The moment-factoring property of 
Gaussian states allow us to write the photocurrent cross- 
correlation at each p 1 as follows, 



C(Pi) 



Co + C n 



dp\K[ n 2 \ Pl ,p)\ 2 \T(p)f 



C p I dp\K { V 2 \ Pl ,p)nT{p)\< 



(5) 



where Co, C n and C p are positive constants, with the for- 
mer constituting a non-image-bearing background pro- 
portional to the product of the incident photon fluxes, 
and the latter two depending on the temporal cross cor- 
relations Ri 2 (t) and R^lir), respectively. The ghost 
image is therefore a linear, space- varying filtered version 
of the object's intensity transmission, whose filter re- 
sponse depends on the spatial cross-correlations between 
Ei and E 2 . Because any pair of phase-insensitive and 
phase-sensitive cross correlations may be associated with 
a classical Gaussian state (i.e., a Gaussian state with a 
proper P-representation), ghost-image formation is in- 
trinsically classical [ill ]. 

Suppose the signal and reference fields have the same 
Gaussian-Schell model phase-insensitive auto-correlation 

K^ m { Pl ,P 2 )RtUh-ti) = 



2P 



e -(|Pi| 2 + |P2l 2 )/io-|P 2 -Pil 2 /2Po e -(*2-t 1 ) 2 /2T 2 



where ao denotes the beam radius, po is the coher- 
ence length (assumed to satisfy the low-coherence con- 
dition po <C ao), and To is the coherence time of the 
P-photons/sec signal (m — S) and reference (to = R) 
fields. Two-mode Gaussian state sources with no phase- 
sensitive cross correlation between the signal and refer- 
ence modes (i.e., Kg R = 0) always have a proper P- 
representation, and hence are classical. A well-known 
example is thermalized laser light, generated by passing 
a continuous- wave laser beam through a rotating ground- 
glass diffuser followed by a 50/50 beam splitter. These 
two beams have maximum \ (Eg(p 1 ,ti)EF>(p 2 ,t 2 )}\, given 
by (J6j) . Taking the phase of this cross correlation to be 
zero we get 

C( Pl ) = C + C„(2P/^) 2 e- 2 l^l 2 / Q « 

x f dpe-l p i-"l 2 / p o e - 2|p|2/a «|T(p)| 2 (7) 

in the near field of the source (i.e., L such that the Fres- 
nel number product Dq = koaopo/2L 3> 1). From ([7]) we 



see that the phase-insensitive cross correlation between 
Ei and E 2 has three important consequences. First, the 
ghost image is space-limited by the average intensity pro- 
file of E 2 , which affirms that the object must be placed 
in the field of view ao of the reference beam 12j . Second, 
the finite size of E\ restricts the useful transverse scan- 
ning range of the pinhole detector to a radius ao. Finally, 
and most importantly, the finite cross-correlation coher- 
ence length po limits the resolution of the image. When 
space-limiting and finite detector area can be neglected, 
the ghost image in Q is proportional to the convolu- 
tion of the object's intensity transmission, |T(p)| 2 , with 
the Gaussian point-spread function e~' p ' / p °. Thus the 
resolution — radius to the e _2 -level in the point-spread 
function — is \[2p§. 

If the object is in the source's far-field (Dq -C 1), the 
image still has the form ([7]), but the beam diameter and 
the coherence length diffract in equal proportion, so that 
ao and po must be replaced with = 2L/kopo and 
Pl = 2L/k Q a Q respectively [l3|. Therefore, the far-field 
field-of-view increases to a^ while the image resolution 
degrades to 2y/2L/k ao. Thus, so long as field-of-view 
is not the limiting factor, it is more desirable to place 
the object in the near field of the source. When the ob- 
ject is far from the source, the source coherence can be 
transferred to the object plane with lenses. 

When Gaussian-state signal and reference fields have 
a non-zero phase-sensitive cross correlation, their joint 
state need not have a proper P-representation, viz., the 
state may be non-classical. Consider signal and ref- 
erence fields with Schcll-modcl phase-insensitive auto- 
correlations given by A*(xi)A(x 2 )g( nS> (x 2 — xi), for x = 
(p, i) , a phase-sensitive cross correlation of similar form, 
A(xi)A(x 2 )g( p } (x 2 — xi), and no phase-insensitive cross 
correlation. Here we shall require |A(x)| < 1, so that 
this function may be regarded as a (possibly complex- 
valued) spatio-temporal attenuation of two homogenous 
and stationary random fields with phase- insensitive auto- 
correlations g^ (x2— xi), and a phase-sensitive cross cor- 
relation g( p,> (x 2 — Xi). Then, the phase-sensitive cross 
correlation spectrum, given by the 3D Fourier transform 
g(p)(f) = JP{ 5 (p)(x)}, must satisfy 



|5 (p) (f)| < ^W(f)(l+gW(f)) 



(8) 



where g( n '(f) = !F{g( n '(x)} is an even phase-insensitive 
auto-correlation spectrum. However, a classical Gaus- 
sian state has its g( p '(f) limited by the more restrictive 
necessary and sufficient condition 



^) (f) | <p(»)(f). 



(9) 



Applying ((9]) to the Gaussian-Schell model auto-cor- 
relations in (|6|) and taking the phase to be zero, we 
find that the maximum classical (Es(Pi,ti)E]i(p 2 ,t 2 )) 
is also equal to The photocurrent correlation C(p 1 ) 
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in the source's near field is found to coincide with the 
thermal light case in ([7]), yielding the same field-of-view 
and resolution. However, phase-sensitive coherence prop- 
agates differently than phase-insensitive coherence [1 51 ] . 
In particular, deep in the far-field — when Dq <^ 1 and 
k al/2L -C 1 — we get, 

C( Pl ) = C + C p (2P/na 2 L ) 2 e- 2 ^ 2 / a i 
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FIG. 2: (Color online) Ghost imaging setup with relay optics. 



where = 2L/kopo and pl = 2L/fcgag are the same as 
the thermal light case. Therefore, (fT0|) has the same res- 
olution and field-of-view as the far- field image from ther- 
mal light, but the image is now inverted, i.e., the ghost 
image is proportional to the convolution of \T{— p)\ 2 with 

The maximum non-classical phase-sensitive cross cor- 
relation satisfies ([8]) with equality. We will restrict our 
attention to two limiting cases in which such a phase- 
sensitive cross correlation is coherence separable, so that 
we may utilize the machinery developed earlier in this 
paper. If g( n \f) 1, the distinction between quan- 
tum and classical sources diminish and the classical 
phase-sensitive results apply. On the other hand, if 
ff (n) (f) « 1, then |3 (p) (f)l w v^W(£), and the source 
exhibits strongly non-classical behavior. Assuming a real 
£T p '(f) and using ©, this approximation holds when 
PT p 2 /a 2 < 1, yielding 



{E s {p 1 ,t 1 )E R {p 2 ,t 2 )) oc 



|Pl| 2 + |P2| 2 )/^-|p 2 -p 1 | 2 /Po p -(*2-tl) 2 /T 2 



(11) 



Equation (fTTj) is a Gaussian-Schell model correlation, so 
C( Pl ) = C + C' p e- 2 ^ 2 ^ 



dpe 



-2|P!-p| 



/Po e -2|p| 



/Q «|T(p)| 2 (12) 



in the source's near field, where C' is a positive con- 
stant. Here the ghost image has field-of-view ag and res- 
olution po, which is a factor-of--\/2 better than the resolu- 
tion from (phase-sensitive or phase-insensitive) classical 
sources with identical source auto-correlations. 

The far-field C(p 1 ) for the non-classical source is given 
by (fTO)) , with p L = 2L/k a , ah = 2\^2L/k po, and C' p 
replacing C p . Hence the far-field resolution for the quan- 
tum source equals that of the classical sources consid- 
ered earlier, but the field of view az,, is %/2-larger. It 
is worth pointing out that these enhancements are due 
to the broadening of the weak spectrum, g( n \{), when 
its square-root is taken. Thus the enhancement factor 
depends on the shape of the cross-correlation spectrum. 



DISCUSSION 

The underlying principle in ghost imaging, whether the 
source is classical or quantum, is the non-zero intensity 
cross-correlation between the signal and reference beams. 
For Gaussian-state sources, this intensity cross correla- 
tion is given in terms of the second-order correlation 
functions, and both phase-sensitive and phase- insensitive 
cross correlations can contribute to the image. Any pair 
of phase-sensitive and phase-insensitive cross-correlation 
functions can be obtained with two classical Gaussian- 
state fields, so long as there are no restrictions on these 
fields' auto-correlation functions. In this respect, the 
ghost image does not contain any quantum signature per 
se. However, if we compare sources that have identi- 
cal auto-correlation functions, we find that non-classical 
fields with low brightness and maximum phase-sensitive 
cross correlation offer a spatial resolution advantage in 
the source's near field and a field-of-view extension in its 
far field. The field-of-view in the near field and the res- 
olution in the far-field are determined by the beam sizes 
at the source, and hence are identical for classical and 
non-classical fields. 

It is worth addressing the physical implications of 
the PT$p\ja\ <C 1 condition, which defines the regime 
in which the aforementioned advantages of non-classical 
light are most prominent. A field that is constrained to a 
radius-ao cross-section and a spatial-bandwidth l/po con- 
sists of approximately ag/pg independent spatial modes, 
and the spectral mode density for a stationary field with 
coherence time Tg is given by 1/T . Hence, PTgpg/ag is 
the source's average photon number per spatio-temporal 
mode, which must be much less than one for (fTTj) to 
be an accurate approximation. However, the source 
has low spatial-coherence, Pq/oq <C 1, so that the to- 
tal number of photons it emits per coherence time can 
be large, PTq ^>> 1, while still satisfying the condition 
PT$p\la\ <C 1. Therefore, the resolution and field- 
of-view improvements achieved with non-classical light 
do not require operation in the biphoton limit, wherein 
each signal-reference photon pair can be time resolved 
by MHz-bandwidth photon-counting detectors. 

Our analysis thus far assumed that the detector plane 
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coincides with the object plane, but a realistic ghost- 
imaging scenario will likely require a separation between 
the two planes, as shown in Fig. O Here the bucket de- 
tector is along an L^-m-long free-space path from the 
object that is not under our control, but the signal- 
arm path may be modified freely. Thus we place a 
focal-length / lens in the signal path at a distance d\- 
m from the object plane and d 2 -Ta from the detector 
plane, where l/d\ + l/cfe = 1//- In addition, because 
the optical path lengths may be different, we introduce a 
post-detection electronic delay (Lr — d% — d 2 )/c to max- 
imize the temporal coherence of the two fields. The 
photocurrent cross-correlation is then C'{p^) = Co + 



j A dp 2 [C n \K^> 2 ,(Pi,P2)\ 2 + C p \K^ 2 ,( Pl ,p 2 )\% where 

Ky^'iPii P2)' f° r m e { n '-P}i are the propagated phase- 
insensitive and phase-sensitive cross correlations whose 
magnitudes are 



k M 



2nL, 



dp'e 



- i 2Tt( 2 P2-p'-|p' 



K[ m 2 '(Mp llP ')T(p') 



with M = —d 2 /d\ being the signal-arm magnifica- 
tion factor. For a sufficiently large bucket detector we 
can utilize Parseval's theorem to show that C'(pi) = 
M 2 C(Mp 1 ), where C(p 1 ) is given by (J5]). Thus, choos- 
ing d\ = d 2 — 2/ yields an inverted version of the object- 
plane ghost image. Image resolution and field of view 
are then determined by the phase-sensitive and phase- 
insensitive coherence properties of the object-plane fields, 
and the placement of the detectors relative to this plane 
only determines the signal-arm optics that are needed to 
obtain this object-plane ghost image. 

So far we have concentrated on the image-bearing 
terms in ([5|). However a background term, Co, is 
also present and it degrades image contrast. To as- 
sess this effect, let us assume that both detectors 
have Gaussian shaped post-detection filters with effec- 
tive integration time Td, and that the temporal part 
of the sources' cross-correlation functions are given by 
([6|) for the classical fields and by (fTTj) for the non- 
classical fields. Then, for classical fields — either phase- 
insensitive or phase-sensitive — the contrast is propor- 
tional to l/yl + T|/4Tq , which is close to unity when 
the fields are narrowband, i.e. when To 3> T<j. How- 
ever for broadband fields, the contrast is proportional to 
the ratio of the detector bandwidth to the source band- 
width, T /Td -C 1, and the degradation is severe. On the 
other hand, for non-classical low-brightness fields — as oc- 
curs in the biphoton limit — the contrast is proportional 
to (al/PTopD/^l + T}/2T%, so that high contrast is 
possible even with broadband fields. 

It is relevant to point out that some features of bipho- 
ton ghost imaging which have been previously associated 
with the entanglement between the two source photons, 
are fundamentally due to the phase-sensitive cross corre- 
lation at the source, rather than the entanglement per se. 



So, when ghost imaging is performed with phase-sensitive 
light, image inversion occurs for both classical and quan- 
tum sources. This is due to the difference in free-space 
propagation of phase-sensitive and phase-insensitive cor- 
relations 15[; it is not necessary for the phase-sensitive 
coherence to be stronger than classical. 

In summary, ghost imaging with Gaussian-state 
sources is due to phase-sensitive and phase-insensitive 
cross correlations, with only the image contrast being af- 
fected by whether the light is classical or quantum. When 
ghost images from classical and quantum sources hav- 
ing identical auto-correlations are compared, the non- 
classical sources with much less than one photon per 
spatio-temporal mode offer resolution enhancement in 
near-field operation and field-of-view enhancement in far 
field operation, in addition to higher contrast in both 
regimes. The large number of spatial modes implies 
that the total photon flux need not be time resolved by 
photon-counting detectors to reap these advantages. Far- 
field spatial resolution and the near-field field-of-view are 
determined by the beam sizes at the source, so they are 
identical for classical and quantum sources. 

This work was supported by the U. S. Army Research 
Office MUM Grant W911NF-05- 1-0197. 
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